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Abstract. We develop rigorous, analytic techniques to study the behaviour of biased 
random walks on combs. This enables us to calculate exactly the spectral dimension 
of random comb ensembles for any bias scenario in the teeth or spine. Two specific 
examples of random comb ensembles are discussed; the random comb with nonzero 
probability of an infinitely long tooth at each vertex on the spine and the random comb 
with a power law distribution of tooth lengths. We also analyze transport properties 
along the spine for these probability measures. 
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1. Introduction 

The behaviour of random walks on random combs is of interest from a number of points 
of view. Condensed matter physicists have studied such structures because they serve 
as a model for diffusion in more complicated fractals and percolation clusters [H [2j [3j 0] . 
In the context of quantum gravity, random combs are a tractable example of a random 
manifold ensemble and understanding their geometric properties can provide insight 
into higher dimensional problems [U [6], [7] . Most of the literature concerns approximate 
analytical techniques and numerical solutions, although there are exact calculations of 
leading order behaviour in some cases [Sj. To this end, it is desirable to have rigorous 
methods for determining the geometric quantities of interest and that is the purpose of 
this paper. 

One such quantity is the dimensionality of the ensemble. On a sufficiently smooth 
manifold all definitions of dimension will agree, but for fractal geometries like random 
combs this is not necessarily true. The spectral dimension is defined to be d s provided 
the ensemble average probability of a random walker being back at the origin at time t, 
takes the asymptotic form t~ da / 2 . This concept of dimension does not in general agree 
with the Hausdorff dimension djj, which is defined when the expectation value of the 
volume enclosed within a geodesic distance R from a marked point scales like R dn as 
R — > oo. 
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We know that for diffusion on regular structures the mean square displacement at 
large times is proportional to t, but for a fractal substrate there is anomalous diffusion 
and the mean square displacement behaves like t 2 / dw , where d w represents the fractal 
dimension of the walk and depends sensitively on the nature of the random structure. 

Biased random walks on combs have also been studied in connection with disordered 
materials, since such a system is a paradigm for diffusion on fractal structures in the 
presence of an applied field [HI [10] . As we discuss later there are several different bias 
regimes. Topological bias, where at every vertex in the comb there is an increased 
probability of moving away from the origin was first studied for a random comb with a 
power-law distribution of tooth lengths in [TT] . Other works have discussed the effects 
of bias away from the origin only in the teeth [12J and only in the spine [13J. The effect 
of going into the teeth can be viewed as creating a waiting time for the walk along the 
spine; the distribution of the waiting time depends on both the bias and the length of 
the teeth and the outcome is the result of subtle interplay between the two. 

In [H] some new, rigorous techniques were developed to study random walks on 
combs. This enabled an exact, but very simple calculation of the spectral dimension of 
random combs. The principal idea is to split both random combs and random walks 
into subsets that give either strictly controllable or exponentially decaying contributions 
to the calculation of physical characteristics. These methods were later reinforced to 
prove that the spectral dimension of generic infinite tree ensembles is 4/3 [T5l [IB] . In 
this paper we use and extend the techniques of [H] to deal with biased walks on combs. 
Some of our results are new; some qualify statements made in the literature; and some 
merely confirm results already derived by other, usually less rigorous, methods. 

The random combs, the bias scenario, some useful generating functions and the 
critical exponents are defined in the next section. In Section 3 we introduce some 
deterministic combs, discuss general properties of the generating functions and establish 
bounds that will be instrumental when studying random ensembles. Section 4 looks at 
regions of bias where the large time behaviour is independent of the comb ensemble or 
simply dependent on the expectation value of the first return generating function in the 
teeth. In Section 5 we compute the spectral dimension in regions of bias where it is 
influenced by the probability measure on the teeth. Two specific cases are considered: 
the random comb with nonzero probability of an infinitely long tooth at each vertex on 
the spine and the random comb with a power law distribution of tooth lengths. Section 
6 examines transport properties along the spine for these same probability measures and 
in the final section we review the main results, compare with the literature and discuss 
their significance. Some exact calculations and proofs omitted from the main text are 
outlined in the appendices. 

2. Definitions 

Wherever possible we use the definitions and notation of [14] ; we repeat them here for 
the reader's convenience but mostly refer back to [H] for proofs and derived properties. 
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2.1. Random combs 

Let Noo denote the nonnegative integers regarded as a graph so that n has the neighbours 
nil except for which only has 1 as a neighbour. Let Nt be the integers 0,1, ... ,£ 
regarded as a graph so that each integer n £ N e has two neighbours nil except for 
and I which only have one neighbour, 1 and i — 1, respectively. A comb C is an infinite 



o 



Figure 1. A comb. 



rooted tree-graph with a special subgraph S called the spine which is isomorphic to 
with the root, which we denote r, at n — 0. At each vertex of S, except the root r, there 
is attached by their endpoint one of the graphs Ng or A^. The linear graphs attached 
to the spine are called the teeth of the comb, see figure 1. We will denote by T n the 
tooth attached to the vertex n on S, and by Ck the comb obtained by removing the links 
(0, 1), . . . , (k — 1, k), the teeth T\, . . . , T& and relabelling the remaining vertices on the 
spine in the obvious way. An arbitrary comb is specified by a list of its teeth {T\, . . .} 
and |Tfc| denotes the length of the tooth. Note that we have excluded the possibility 
of a tooth of zero length. This is for technical convenience in what follows and can be 
relaxed [T7] . 

In this paper we are interested in random combs for which the length I of each 
tooth is identically and independently distributed with probability This induces a 
probability measure \i on the positive integers and expectation values with respect to 
this measure will be denoted (•) ft . In particular we will consider the two measures 



i4 



vf = i£> a>1 - W 
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However, the results proved for apply to any measure with the same behaviour at 
large I and we note in passing that the methods used here will work for any distribution 
that is reasonably smooth, for example the exponential distribution. The measure fi B 
has been discussed quite extensively in the literature but \i A has not. 



2.2. Biased random walks 

We regard time as integer valued and consider a walker who makes one step on the 
graph for each unit time interval. If the walker is at the root or at the end-point of a 
tooth then she leaves with probability 1. If at any other vertex the probabilities are 
parametrized by two numbers £\ and £2 as shown in figure 2a and the allowed range of 
these parameters is shown in figure 2b. For walks in the teeth there is bias away from or 
towards the spine depending on whether £ 2 is positive or negative; similarly a walk on 
the spine is biased away from or towards the root depending on whether e\ is positive 
or negative. When there is no bias we say that the walk is 'critical'; the fully critical 
case £1 = 62 = was covered in [T3]. The notation 

b- = 1 - £\ - £ 2 , 
b+ = l + ei- £ 2 , 

b T = l + 2e 2 , (2) 

will be used where applicable since these combinations appear often in our analysis. We 
denote by B, B', Bi, B 2 etc constants which depend on ei and £2 and may vary from line 
to line but are positive and finite on the relevant range; other constants will be denoted 
c, d etc. 

The generating function for the probability pc{t) that the walker on C is back at 
the root at time t having left it at t — is defined by 

00 

Qc(x) = J2( 1 -x) t/2 Pc(t). (3) 
t=o 

Letting uo be a walk on C starting at r, u(t) the vertex where the walker is to be found 
at time t, and the probability for the walker to step from uj(t) to u(t +1), we have 

H-i 

Qc{x)= Up^ty (4) 

w.r— *r t=0 

A similar relation gives the generating function for probabilities for first return to the 
root, Pc(x), except that the trivial walk of duration is excluded. The two functions 
are related by 

and it is straightforward to show that Pc(x) satisfies the recurrence relation 

Pc{x) = (6) 
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Figure 2. Bias parameterisation. 



Note that Pc(x) and Qc( x ) depend upon ei and t2\ to avoid clutter we will normally 
suppress this dependence but if necessary it will appear as superscripts. It is important 
for what follows that Qc is a convex function of Pc which is itself a convex function of 
P Tl ,...P Tk ,P Ck for any k > 0. 

For an ensemble of combs, we will denote the expectation values of the generating 
functions for return and first return probabilities as 

Q(x) = (Qcix)), 

P(x)=(P c (x)) lt . (7) 

We will say that g(x) ~ f(x) if there exist positive constants c, c', a, a' and x such 
that 



cf(x)exp {-a{\og\f{x)\) 1 / a )<g{x)<d f{x)\\ogf{x)\"' 



(8) 



for < x < xq. The tactic of this paper is to prove bounds of this form for the generating 
functions; in almost all cases our results are in fact a little stronger having a = a' = 
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when we will say that g(x) ~ f{x). 

The random walk on C is recurrent if -Pc(O) = 1 in which case we define the 
exponent (3 through 

1 - P c (x) ~ x p . (9) 
If (3 is an integer then we expect logarithmic corrections and define j3 if 

l-Pc(x) ^x^hgx]-^. (10) 
It follows that Qc(%) diverges as x — > and we define a by 

Qc(x)~ x - a , (11) 
and if a is an integer, a when 

Q c (x) ^ x~ a \ log xf. (12) 

If -Pc(O) < 1 then the random walk is non- recurrent, or transient, and Qc{%) is finite 
as x —>■ 0. Then if, as x — > 0, the first fc — 1 derivatives of Qc(a;) are finite but the kth 
derivative diverges we define the exponent by 

Q%\x) ~x~ af % (13) 

and if is an integer, when 

Q { c\x) ^ x- ak \\ogx\ &k . (14) 

In considering the ensemble of combs /i, we define all these exponents in exactly 
the same way simply replacing Pc(%) with (Pc{x))^ and so on. Note that for a single 
recurrent comb f3 = a but in an ensemble this is no longer necessarily the case; applying 
Jensen's inequality to (JSj) we see that j3 < a. 

If QW(i) ~ x~ ak then it is straightforward to show that 
a- 1 

R k (X) = J2t k (PcW^-X-**. (15) 
t=o 

It follows that if the sequence decays uniformly at large t, which we do not prove, then 
it falls off as t OLk ~ 1 ~ k . Thus we define d s = 2(1 + k — af.). Similarly if Q(x) ~ | logx|° 
then R(X) ~ | log Aj" and, again assuming uniformity, pit) falls off as t" 1 \ logtl 0-1 . 

2. 3. Two-point functions 

Let Pc(t; n) denote the probability that the walker on C, having left r at t — and not 
subsequently returned there, is at point n on the spine at time t. The corresponding 
generating function, which we will call the two-point function, is defined by 

oo 

G c (x-n) = Y,^-x) t/2 Vc(t;n). (16) 
t=o 

Letting u be a walk on C starting at r and ending at n without returning to r we have 

|oj|-1 

G c (x;n)= £ {1 - x)M f[ p u{t) . (17) 
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Following the discussion in section 2.2 of [14] this leads us to the representation 

n-i , 

= h n _ ^ II r Pc >W- 



T/ie i/ecrf kernel 

Let Kc{t; n, £) denote the probability that the walker on C, having left r at t = 0, is at 
point £ in tooth T n at time t. Kc{t; n, £) satisfies the diffusion equation on C so we call 
it the heat kernel. The probability that the walker has travelled a distance n along the 
spine at time t is given by 

oo 

K c (t;n) = Y,Kc(t;n,e), (19) 
and has generating function 

oo 

H c (x; n) = ^(1 - xf 2 K c (t; n). (20) 



*=o 



Hc(x; n) can be written as 



Hc(x-,n)= Gc ^f D ]Tnl (x), (21) 
1 — r c {x) 



where 

DJx) = 1 4- 

3 

and we define 



h 

Di{x) = l + -Z-^G Nl (x;k), (22) 



H(x;n) = (H a (x;n)) IM . (23) 
Note that, because Ko(t;n) is a probability, 



5>(x;n) = i±^Hf (24) 



n=0 

The exponent c4 is defined through the moments in n 

oo 

^n k H(x;n) ^ x~ l - d \ (25) 

n=0 

and in the case dk = the exponent c4 is defined when 

oo 

^n k H(x;n) ^ x- 1 \ log x\ dk . (26) 

ra=0 

If €i > one can show that on any comb (w) w: ia,|=f is a non-decreasing sequence and 
thus that there is some constant T such that for T > T 

^ (t£t\) < (W^K + ^M^+i),, < ^ Tdl ' (27) 
c | logT| Jl < ((7i) wsM =T + (n) w:M= r + i)^ < c|logT| d \ dt = 0. 
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If e\ < (for which we always have d\ = 0) then we have only the weaker result that 
for T > T 

(rp \ 1+di T 
n-n) <E« n WN>M <5Tl+ *- ^ 
1 S 1 ' t=0 

3. Basic properties 

5.1. Results for simple regular combs 

The relation ([6]) can be used to compute the generating functions for a number of simple 
regular graphs which will be important in our subsequent analysis |14j . 

(i) An infinitely long tooth, iV^: 

1 — x^ if e 2 = 0; 

P 00 (x) = { i_ 2 |e 2 | x , t o , a , . . (29) 

i (1 — 2e 2 ) + U(x ) otherwise. 



bx 4 e 2 



(ii) A tooth of length £, Nf 



where 



1 -l- XY 1 ^ * 

P e (x) = f w (x) 1 + xy _ t (30) 

brjl-P^x)) Y _ 2-b T P 00 {x) 
2-6 T (l + P 00 (x))' frrPooCx) ' 1 J 

(iii) The comb Jj given by {T^ = N\,\/k} has all teeth of length 1, and 

{1 - B lX ^ + 0(x) if ei = 0; 

'- r N_ ft +0(ja) otherwise . (32) 
o+ |ei| 

Note that (j is non-recurrent if ej > 0. It is also convenient to define 1% to be 
{T 1 = N t ,C 1 = ^}. 

(iv) The comb * given by {T^ = Noo, \fk} has all teeth of length oo and is non-recurrent 
for e 2 > 0, 



Otherwise 



PAx) 



< \ N __g2 if e 2 = 0, ei ^0; 

1 + ei |ei| 

1 — e 2 — jei| 5 3 2 (34) 

; a^i — r + C(x ) if e 2 < 0, ei 7^ 0; v ; 

o+ |ei| 

k 1 - B 4 x^ + 0(x) if e 2 < 0, ei = 0. 
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(v) The comb \>l given by {T^ = Np, \/k} has all teeth of length i and 



P M (x) 



^ - + 1 + h\) x + 0{xH 2 ) if e 2 = 0, e x ± 0; 

1 + ei |ei| 

l ~ e2 ~ l ei l _ + 0{xY- i ) if e 2 < 0, Cl ^ 0; (35) 

o+ |eie 2 | 

1 - — ^ + 0{x*Y- 1 ) if e 2 < 0, ei = 0; 

M 

where, as x — > 0, 

i — ^|e 2 | 

When e 2 > let I = [\ loga;|/logFj, where [z\ denotes the integer below z. For 
£ > 2£ the teeth are long enough that -Pw( x ) behaves like (1531 . For I < £ < 2£, 
P{,e{x) is non- recurrent with the leading power of x being fractional. For I < I 

^M-x^ + 0(x) ife^O; 
H<i{x) = { b ■ 6162 - ( 37 ) 



+ 0{x*Y-* t ,xY- t ) if ei = 0, 
where the notation 0(a,b) means 0(max (a, &)). 

3.2. General properties of the generating functions 

The generating functions for any comb satisfy three simple properties which can be 
derived from (EJ): 

(i) Monotonicity The value of Pc(%) decreases monotonically if the length of a tooth 
is increased. 

(ii) Rearrangement If the comb C is created from C by swapping the adjacent teeth 
T k and T k+1 then P c >(x) > P c {x) if |T fc+ i| < \T k \. 

(iii) Inheritance If walks on C k or T k are non-recurrent for finite k then walks on C are 
non-recurrent. 

The proof of the first two follows that given in [2] for the special case e 2 = e\ = 0. The 
third can be shown by assuming that either -Pci(O) < 1 or -Pr\(0) < 1; it then follows 
immediately from (jfJJ) that -Pc(O) < 1 and the result follows by induction. 

3.3. Useful elementary bounds 

By monotonicity Gc(x; n) is always bounded above by G$(x; n) from which we get 

G c (x;n) < — exp(-nA ei ' e2 (x)), (38) 

0-4- 
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where 



A ei ' e2 (x) 



2 + e 2 

26! 



X 2 



log 



2 + e 2 
l-e 2 



if e x > 0, 

if ex = 0, 
if e l < 0. 



(39) 



Now let Pq(x) denote the contribution to Pc(x) from walks that reach beyond n = N 
on the spine. It is straightforward to show using the arguments of section 2.5 of [bi 
that 



P^ N \x) < ~b-G^ 2 (x;N)Gc ei ' e2 (x;N). 



1. 



Combining this with 



we obtain the useful bound 
3b. 



Pc N \x) < ^exp(-iV(A ei ' e2 (x) + A- ei ' e2 (x))). 



(40) 



(41) 



Now consider the ensemble y! of combs C for which: = iVi, /c = 1..K — 1, 
Tft- = A^; at k > K teeth are short, — N%, with probability 1 — p or long, = A^, 
with probability p; and the nth tooth is short, T n = A^. Then using the representation 
(JT8J) Gc(s; n) can be bounded above by noting that if T k+1 = N e then Pc fc < P^, 



otherwise Pc fc < Pjj. This gives 



Gc(x,n)< ^J'p^ 



.k+K 



(42) 



and hence 
(G c (x,n)) , 



< 



n-K-l 

£ 

fe=0 

3 A 



n-K-l 



P 



n-K-l-ki 



p) k G c (x,n) 



P^{x) K P n {x) 
(I - xY' 2 



((l-p)P $ (x)+ P P ei (x)) 



n-K-l 



(43) 



4. Results independent of the comb ensemble ji 



In this section we show that in some regions of ei >2 the behaviour at large time is 
essentially independent of the comb ensemble, or else simply dependent upon (Py(rr)) . 
The leading, and where different, the leading non-analytic, behaviour of (Pr(x)) as 
x — > for the measures studied here is given in table [TJ The results for \i A are trivial, 
as are those for any measure when e 2 < 0, while the case y B and e 2 = can be derived 
using the techniques in [H]. The calculation for y B and e 2 > is somewhat subtle and 



is included in Appendix A 
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Table 1. Leading and leading non-analytic behaviour of 1 — (Pt) u in various cases. 



ensemble 


e 2 < 


e 2 = 




e 2 > 








Bx 


Bx^ 




B + B' 


t 

X 




[i B ,a<2 


Bx 


Bx a/2 




5(1 log 


x\ 


a-l\-l 


[i B ,a = 2k 


Bx 


Bx + . . 


.B'x k \\ogx\ 


5(1 log 


x\ 


a-l\-l 


fi B ,a>2,«^ 2k 


Bx 


Bx + .. 


. B'x a / 2 


B(\ log 


x\ 


a-l\-l 



4-1. d s when e 2 < 

First we show that for any comb ensemble 

if e 1 < and e 2 < 0; 

1 if ex — and e 2 < 0. 
By monotonicity we have that for any comb C 

P*(x) < Pc{x) < Pt(x). 
Taking expectation values and using (1321) and 

1 - B x x^ + 0{x) 
l-e 2 



P(x) = {P c {x)) L 



X- r + U(X 

N 



it follows that for e 2 < 
if ei = 0, 
otherwise. 



Similarly 



and so 



Q*{x) < Qc{x) < Q$(x) 



^+0(1) 

J2 



B 2 \e 1 \ 



0(1) 



x 



if ei = 0, 
if ei < 0, 



(44) 



(45) 



(46) 



(47) 



(48) 



Q(x) = {Qc{x))^ - 

and (T44]) follows. 
4-2. d s when e\ > 

When ei > all combs are non-recurrent and so we must examine the derivatives of 
Q(x). Differentiating (j5J) and flBj) gives 

Qg ) (x) = g c (x) 2 p( 1) (x), (49) 



^(z) 



X 



Pc(x) 2 
1 - x)6_ 



b T P^(x) 



b + P$(x 



(50) 
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By monotonicity (l5Tf|) can be bounded above and below by replacing Pq with and 
Pjl respectively. Taking the expectation value and using translation invariance to note 
that (Pc)n = (Pc-l) u shows that, if \Prp (x)\ diverges as x — > 0, then 



~ B (P$\x)) + 0(v^ (p?\x)) > 1)- (51) 

As can be seen from tabled! in some cases (Pt(x)) is analytic, or only higher derivatives 
diverge. For the measures considered here it can be shown that if (Pt{x))^ is analytic 
at x = then so is Q(x). If on the other hand (Pt(x)) is not analytic but the fc'th 
derivative diverges then 

Q^(x) = B (4 fe) (z)) + +0(v^ (P?\x)) , 1). (52) 
The proof is a straightforward but tedious generalization of ( )49l) and ( 1501) and is relegated 



to Appendix B If a derivative of Q(x) diverges then d s can be read off using ({14")) and 
( I52p . Otherwise if all finite order derivatives are finite then pc{t) decays at large t faster 
than any power and d s is not defined. 

4-3. dk when e 2 < or e\ < 

We show that for any comb ensemble 

{0 if ei < 0, 

k/2 if e x = and e 2 < 0, (53) 
k if ei > and e 2 < 0. 

It is trivial to show that 

1 < D e < ^- e 2 < 0, (54) 

l e 2| 

and then by monotonicity we get 

< H(x\n) < - — | — . (55) 



l-P*(x) ~ v ' J ~ \e 2 \l-P t (x 

Combining this with ( |32i) and ( |34|) yields the results for e 2 < 0. 
To deal with e\ < and e 2 > note that monotonicity gives 

f lT p X ( \ ) G m (x;n)<H{x;n)<( D ™ { *\ ) G t {x;n). (56) 
1-P c {x)/^ \1-P (x)/^ 

Using the lower bound and (fT8|) . (1241) and (1331) we get after summing over n 

3 ''- <±H(x;n)< 2 -. (57) 



1 - P c (x) I ^ 6+ ^462 + ef - ei - 2e 2 ^ a; 
Inserting this into the upper bound of (|56[) gives 



H(x; n) < l^ V^+^-ei- ^ ^ ^ (5g) 
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It is a trivial consequence of (1241) that 

oo 

^n k H{x-n) > -, k>0, (59) 

ra=0 

and the results then follow by using ( 15B1) . 



5. The spectral dimension when e 2 > and ei < 

Here and in some of the sections to follow we will need to sum over the location of the 
first long tooth to determine the spectral dimension. Most generally we call a tooth 
long when it has length > £ and short when it has length < £. Consider combs for 
which the first L — 1 teeth are short but the Lth tooth is long; the probability for this 
is p{l — p) L_1 , where p is the probability of a tooth being long. Denoting by £L a comb 
having the first long tooth at vertex L gives 

oo 

Q{x)=Y t {Q IL {x)) ll p{l-p) L -\ (60) 

L=l 

Qil(x) is bounded above by the comb in which all teeth at n > L + 1 are short, and 
below by the comb in which all teeth at n > L + 1 are infinite, 

Q{T n <L=N tl ,V<i;T n > L =N^}{x) < QtL{x) < Q{T njtL =Nv,T L =N e }(x). (61) 

5.1. ji A - Infinite teeth at random locations 

5.1.1. 62 = 0, ei < We first show that the exponent /3 = \ - so it is unchanged from 
the comb *. This result follows from the inequalities 

l-t-r—r + 0(x) <P(x) < 1-pB'x? + 0(x). (62) 
l e i| 

The lower bound is obtained by applying Jensen's inequality to (J6j). To get the upper 
bound we average over the first tooth and then by monotonicity we obtain 

P(x)<pP ei (x) + (l-p)Pz(x), (63) 

with £ = oo and using (jBJ) and (1521) gives the bound required. 
The spectral dimension is given by 

f 1 ifp>2|e 1 |(l + | ei |)- 1 , 

^ s ~ 1 — ^7 ^r- otherwise. 

[Mm) 

This result follows from estimating the sum in ( 1601) using the bounds in ( 16T1) with £ = oo 



and short teeth being N\. Pc(%) for these bounding combs is computed in Appendix C 
and using (1C.3I) we get upper and lower bounds on Qool(^) of the form 

1 (65) 



Bx + B'x* (0) 



L ■ 
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5.1.2. e 2 > 0, e\ < The probability that C is non-recurrent is at least p, the 
probability that T\ = N^, and hence 

P(0) < 1. (66) 

In fact it follows from the lemma of Appendix B that p( k \x) is finite for all finite k so 



the exponent (3 is undefined. 

The spectral dimension is given by 



= 2bg(l-rt (67) 

i og f afa ' 



To show this we start by estimating Q(x) in exactly the same way as in 15.1.11 except 
that the behaviour of the limiting combs is now given by (IC.5I) so that there are upper 
and lower bounds on Qool(%) of the form 

' • (68) 



|ei| — £2 / 

When p < 1 — b + /b^. this sum diverges at x = and it is then straightforward 
to obtain fIBTj) . For larger p the sum is convergent at x = so we next examine 
Q^(x) = (QqPq^J . Note that —P^ > then letting Z be a very large integer 
and using Holder's inequality 

6_ <Qc(*) 2 ) M < (x) < (Qc{x)™' z )f* (-PW(,) 2Z+1 )^ . (69) 



By the lemma of Appendix B the second factor in the upper bound is finite as x — ► 
so we need an estimate of (Qc)u- This i s provided by fl68l) modified by squaring the 
denominator; when p < 1 — (b + /bJ) 2 this sum diverges at x = and once again we 
obtain (1671) . For still larger p both Q and are finite at x = and we examine 
the second and higher derivatives. This uses 

PUD , (-l) fe -Po ) > fc_fc+~ 1 / 32fe_1 ) Hdlder's 
inequality and the lemma; the term with the highest power of Qc dominates and the 
result is always fl67|) . 

5.1.3. e 2 > 0, ei = By the same argument as in l5.1.2l we find P(0) < 1, so (3 is again 
undefined. An upper bound on may be obtained as in 15.1.11 using ( 1C.9I) to get 

QooL^) < (L + (l-e 2 )/4e 2 ) (70) 

which means the upper bound of (160]) is finite. A proof that all derivatives of Q(x) are 



finite is given in Appendix B.2, so pc(f) decays faster than any power at large t. 



\ Strictly speaking when 1 — (b + /b-) k < p < 1 — (6 + /6_) fe+1 / z the upper bounds diverge so our proof 
does not work for these arbitrarily small intervals. 
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5.2. \i - Teeth of random length 

In this subsection we are concerned with random combs that have a distribution of tooth 
lengths. The general strategy for determining quantities of interest is to identify teeth 
that are long enough to affect the critical behaviour of the biased random walk and 
consider the probability with which they occur. It will be useful to define the function 

wc s\ i log^T — C(a — 1) log | loga;| | , 
A (<^ V, = I p J , (71) 

which will be used to denote a tooth length, and the function 



M^fr^c^n, (72) 

which is the probability that a tooth has length greater than I — 1. 



5.2.1. e 2 = 0, ex < We first show that 
f - if a < 2, 

P= 2 (73) 
I 1 otherwise. 

The proof follows the lines described in section 15.1.11 with a slight modification for 
the upper bound on P(x). Note that, from (150]) . teeth of length £ > [x~^\ have 
Pt{x) < 1 — Bx2. We then proceed as in (163]) but with £ = [x~i\ + 1. 

The exponent f3 is non-trivial if a < 2 but, as we now show, d s = for all a > 1 so 
mean field theory does not apply when a < 2. This result follows from the inequalities 

, B \^_ < Q{x) < -• (74) 
xllogxl"- 1 x 

The upper bound is a consequence of Q(x) < Q$(x). To obtain the lower bound consider 
the combs for which at least the first iV teeth are all shorter than £q. Then using 
monotonicity and ( )4Ti) 



0(x) > (ItaM! (75 ) 

Vl ; - l-P Mo (x) + 0(exp(-N(A^(x)+A-^(x)))Y 1 ; 

Setting E = A(l, (a - l)" 1 , 0), N = [2(A ei ' e2 + A -61 ' 62 )^ 1 ] logx|J + 1 and using ([35]) the 
result follows for small enough x. 

5.2.2. e 2 > 0, ex < The exponent (3 — but there are computable logarithmic 

corrections and we find that 

B . . B' 

— r— 7 < P{x) < 1 - 1 r 

logx| a 1 | log a? | 



The lower bound follows from applying Jensen's inequality to ([6]). For the upper bound 
note that teeth of length £ > A(l, 1,0) have Pt < B. Again proceed as in (163]) with 
£ = A(l,l,0) + 1. 
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Table 2. (Dg) in various cases. 



ensemble 


e 2 < 


e 2 = 


e 2 > 










B + 0(x) 


Bx-i +0(1) 


Bx" 1 -+■ 


0(1 


) 




fi B ,a>2 


B + 0(x) 


B + 0(x) 


B(x\ log 




!)-!_ 


1-0(1) 


fi B ,a<2 


B + 0(x) 


Bx a/2-l +0(1) 


£(x log 


;x Q " 




1-0(1) 



The spectral dimension is d s = showing again that mean field theory does not 
apply. This follows from the inequalities 

B'exp(-B"\\ogx\ 1/a ) B , . 
— ; <Q <—, (77) 

x x 

for small enough x. The upper bound is a consequence of Q(x) < Q${x) and the lower 
bound follows from fSJ by setting £ = A(l, 1/a, 0), JV = L2(A £1 ' e2 +A" ei ' e2 ) _1 | logx|J +1 
and using (1371) . 

5.2.3. €2 > 0, E\ = The exponent /3 = 0, but there are logarithmic corrections which 
follow from the inequalities 

The lower bound comes from applying Jensen's inequality to the recurrence relation (JH]). 
The upper bound is obtained by requiring unitarity of the heat kernel and its proof is 



relegated to Appendix D 



The spectral dimension and logarithmic exponent are given by 
4 = 2, 

a = a — 1, (79) 

which shows that mean field theory does not apply. This result follows from 

B' (logxl^ 1 < Q(x) < Bllogxl"- 1 (80) 

for small enough x which is obtained by a modified version of the argument in 15.1.11 
First let £ = A(l, 1, (), so that 

P in (x) = l-- ^ 77 -^ + O f ' 



l ogx |C(a-l) \J log 2 1 ^t^ 1 ) 

To obtain (IHDl we use (!60|) and (!6T!) with p = p > (£o), £ = £o and for the lower bound 
set T n< i = Ne . Then using the bounds in ( 1C.8I) with ( = 1 and (1C.7jl with ( = 2 and 
estimating the sums gives the result. 

6. Heat Kernel when t\ > 0, e 2 > 

These calculations require (D^) in the various cases which are tabulated in table [2] for 
convenience. 
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6.1. fi A - Infinite teeth at random locations 
We show that 

f if e 2 > and d > 0, 

4 = < ~ 82 

[ k/2 if e 2 = and e x > 0. 

These results follow from fl85l) . (IBS"]) and flH7|) below. 

Noting that for ei > all combs have 1 — BZ 1 — < Pc(%) < 1 — -B+ 1 and using 
monotonicity gives 

B- (D^ix))^ G.(s;n) < H(x;n) < (D^x))^ ( f^^ ) , (83) 

<B+(D vrn \{x)) it {G cl {x;n)) li , (84) 

where C is constructed from C by forcing T n = N\. If e 2 > then using fj43l) with 
= 0, £ = oo gives the upper bound 

H(x;n) < — exp(-B'n). (85) 

If e 2 = then exactly the same calculation gives 

B i 

H(x; n) < — exp(— B'x^n) (86) 

and evaluating the left hand side of (|83|) gives a lower bound of the same form. If e 2 > 
and ei = it is necessary to sum over the location of the first infinite tooth. Using 
(IC.9I) . fj43|) and introducing C as in fl83l gives 

H(x;n) < — exp(-B'n). (87) 

x 

6.2. fi B - Teeth of random length 
We show that 

if e 2 > and e x > 0; 

d k = { ka/2 if e 2 = 0, e x > and a < 2; (88) 
A; if e 2 = 0, ei > and a > 2. 

These results follow from 



— -exp(-B'n/\ logx| a_1 ) if e 2 > and ei > 0; 



H(x: n) < < 



x\ logx| 

— — exp(-B'nx a/2 ) if e 2 = 0, e x > and a < 2; (l>9) 



x l-a/2 

k Bexp(-nB'x) if e 2 = 0, ei > and a > 2, 

when x is small enough and lower bounds of the same form. 

The upper bounds are obtained by proceeding as in subsection 16.11 for e% > and 
e 2 > setting £ = A(l, 1, 0) + 1 and for e x > and e 2 = setting £ = [a^ 5 J + 1. For 
ei = and e 2 > we start with the upper bound of (I83I) ; let t\ = A(l, 1, 2), p x = p > (£i) 
and £ 2 = A(2, 1, 0), p 2 = p > (£ 2 ). The latter shall be called long teeth and we denote by 
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(^2-K"tO the comb with a single long tooth at vertex K. We now sum over the location 
of the first long tooth using ffl8|) . fj43|) and (1C.8|) and taking account of the fact that the 
first long tooth may be before or after the nth tooth 

n-l _ \K-X K ~ l 



P2(l-Pl) K -' 

OO M x 

\^ P2(l - Pi) 



K=1 - - fcA-BK-, m=0 

X ((1 - Pl )Ptf + (Pi - P2)^ lS + P2^ 2 fl) n " X_1 

\/sT-l 



In the first sum we use the value given in Table [2] for {D\T n \{x)) ■ In the second sum 
(0(i 2 ~ |T n |)D| Tn |(a:)> = B(x\ logx] 2 ^ 1 ))- 1 + 0(1) for \T n \ < £ 2 and the result follows. 
To obtain the lower bounds when e\ > we note that 

H(x;n)>B„(D lTnl (x)G c (x;n))^ 

= B_(D lTnl (x)) ii (G c (x;n))- (91) 
where the measure fx is defined by 

~p e = fie, for teeth T^, k ^ n, 

= , for tooth T n . (92) 

Using the decomposition (1T81) and Jensen's inequality 

(Gc(x;n))-> 3(1 ~ + x) " /2 exp(-^ n ), (93) 

where 

Sn = £(f -JW*)) +f <l-PT fc+1 (x)>_. (94) 

k=0 ^ + / IM + 

Now applying Jensen's inequality with the measure /Z to (jSJ) shows that the lower bounds 
satisfy a recursion formula of exactly the same form as discussed in Appendix C So 
from ( 1C.1I) we find that 

S n < n (*± - P(x) + ^ (1 - Fr(x)) t 

V P(x)(l-A(x)) 

^(xJfc-^PCx) - A(x)P(x))/(P(x) - P{x)) - 1' 

;i -rc)6_ 



where 



P(x) 
P(x) 



3-6 T (P T (x)) M -6 + P(x) 

(1 - x)6_ 
3-6 T (P T (x))_-6 + P(x) 
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For ei > it is straightforward to check that A(x) > c > 1 and that the sum in (1951) is 
bounded above by an n independent constant. Lower bounds of the form of fl89|) then 
follow by inserting the appropriate (Pt)^ in (EE) and ( 1951 . 
When ei = 

H(x;n)>(D ]Tn] (x)) ( ^ G ' { ^\ ) > ( 97 ) 



where C" is constructed from C by setting T fc > n = N^. Choosing £ = A(l,l,2) and 
using (fT8|) and (1C7|) gives 

H(x;n) > (D lTnl (x)) (1 - P> (4)r 1 3(l -x)-/ 2 ^ 2 ^ f " 

1 -H l 3 '/ i" n -l 

> <%r„|(*)>^ (1 - P>(4))- 1 3(1 - xY^P^xf 

x -i * , exp(^- 3 ^-^^l (98) 

for n > 4 which gives the result. 



7. Results and discussion 

Figure [3] outlines the results that we have computed for \i A . These are new and show 
that the most interesting regime is actually when the bias along the spine is towards 
the origin, a circumstance which has not been studied much in the literature. When 
ei > and 62 > the walker disappears rapidly, never to return, and p(t) decays faster 
than any power. When €\ < the bias along the spine is keeping the walker close 
to the origin but if there are any infinite teeth present the walker can spend a lot of 
time in the teeth; the conflict between these effects leads to a non-trivial d s . The fact 
that d\ = whenever t 2 > shows that the walker never gets far down the spine; if 
she disappears then it is up a tooth that she is lost. The Hausdorff dimension for \i A 
is dn = 2, regardless of bias and so we have here several examples of violation of the 
bound 2dn/(i + dn) < d s < dn, which applies for unbiased diffusion [19] . 

Figure H] shows our results for jj, B as well as the results for the unbiased case studied 
in [13] . This length distribution has been studied quite extensively in the literature but 
usually under the assumption that e 2 > 0. As can be seen the interesting behaviour 
displayed by \i A when E\ < does not occur here - essentially because very long teeth 
are not common enough. We believe that with more work a when e\ < 0, t 2 > can be 
found using our methods, but as this will not give further physical insight we leave the 
calculations to elsewhere [17] ■ The case e± > 0, e 2 > (often called topological bias) was 
originally studied using mean field theory, which gave the mean square displacement 

(n 2 (t)} ~ (logt) 2 ^, (99) 

and this is in fact correct since the walker spends much of the time in the teeth. However 
the claim in [3] that (|99|) holds for e 2 > regardless of e\ is false. The mean field method 
gives the correct result when ei = only because the walk on the spine is ignored, which 
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amounts to using Pt(x) for Pc(x) in (|2"Tj) and naively applying Jensen's inequality. The 
case €\ > 0, e 2 = was studied by Pottier [13] who computed the leading contribution 
exactly, but without complete control over the sub-leading terms; she also calculated 
the leading behaviour (n 2 ) — (n) 2 which we have not. Of course our results for d s and 
d\ agree with hers. The Hausdorff dimension for fi B is dn = 3 — a when a < 2 and 
ds — 1 when a > 2 and so again we see that, as expected, a biasing field intensifies 
the difference between the purely geometric definition of dimension and that which is 
related to particle propagation. 

The results for e 2 < are intuitively obvious and, as we have proved, apply for any 
model with identically and independently distributed tooth lengths. The walker never 
gets far into the tooth and therefore combs have long time behaviour characteristic of 
the spine alone. 

This paper has given a comprehensive treatment of biased random walks on combs 
using rigorous techniques - namely recursion relations for generating functions combined 
with unitarity and monotonicity arguments. It serves to put in context many previous 
results as well as present new ones. In the unbiased case [H] and in some bias regimes 
mean field theory is sufficient to compute the leading order behaviour because the walker 
either does not reach the ends of the longest teeth or does not travel far enough down 
the spine for variations from average to be important. But, as is illustrated in many 
examples here, a full treatment is needed when such fluctuations cannot be ignored. 
Finally, while the results are of interest in themselves, an important point of the paper 
was to demonstrate that rigorous analytic methods can be used to treat biased diffusion 
on random geometric structures and it is to be hoped that these tools can be extended 
to higher dimensional problems. 
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Appendix A. Calculation of {Pt{x)) B for e 2 > 
First we rewrite (!30|) as 



Pi(x) = Poo(x)Y - (Y - l)Poo(x)X 



1 



1 



(A.l) 



x -i + Y -t> 



so that 




(A.2) 
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/ 


£2 


\ 




/ 

/ 

/ 

/ 

yd s = 21og(l -p) IT 1 
X 4 = 


4 n.d. 

4 = o 


N 

N 

s 

\ 

\ 

N 

\ 

N 

\ 

4 n.d. 

4 = 




4 = 1 p > p* 
4 = log(i -p)Q _1 p < p* 
4 = 


d - - 

^7 — k 

a k — 4 
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dk - 2 


ei 


d s = 

4 = 

/ 

/ 

/ 

/ 

/ 
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4 = | 


4 n.d. \ 

4 = k 
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\ 

\ 









Figure 3. Results for /r 4 where £1 = log ^"j^jl^ j and p* = 2|ei|(l + |e x | - e 2 ) _1 . 
The logarithmic exponents a and dfc are always zero for \i A . 
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/ 
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4 = o 
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4 = 

4 = ^(a — 1) 


N 

4 = 2 

a = —a 

4 = 

4 = k(a — 1) 




4 = 
4 = 
4 = 


4 = f if a < 2 
4 = | if a > 2 
4 = if a < 2 
4 = | if a > 2 


4 = 2 + a 

a*; = 1 if a = 2A; 
cifc = if a 7^ 2A; 
4 = f if a < 2 
dk = k if a > 2 
4 = 
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4 = 

/ 

/ 

/ 

/ 

/ 

/ 


4 = 1 

4 = | 


4 n.d. \ 

4 = k 

\ 

\ 

\ 

\ 

\ 
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Figure 4. Results for [i B . When e 2 < the logarithmic exponents a and dk are 
always zero. 
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with 

cj- a 



' v 1=1 



Since for e 2 > 0, Y > 1 we let log F = p and write 



l CT ii£g£ij 



5 " E + e-P« + E X "i a +e-^' (A ' 4) 



£=1 ^^f^J+i 



where a is an arbitrary constant < 1. This is bounded above by taking i in the 
exponential to be its value at the top of each sum to give 

5 < - + rr ^ T - (A.5) 
x CT | log a; | 

Noting that as x — > 0, X -1 — >■ Bx we get a lower bound on (Pg(x)) B of 

A 

logx| 

for small enough x. An equivalent upper bound is calculated in the same manner by 
ignoring the first term in ( 1A.4I) and setting a — 1, which leads to the result quoted in 
tabled! A similar procedure leads to bounds of the form B/x\ logx| a on (^P^\x)^ , 
which we also need, at small enough x. 

Appendix B. Proof of results for non-recurrent regime 

First we define a structure of ordered lists of ordered integers. Let S denote an ordered 
list of hs integers 

([n 1 ,n 2 ,...n hs ], n x > n 2 > . . . > n hs > 1, h s > 1, 
~\[], h s = 0. 

Define 

hs 

, > Hi, hs > 1, 
\S\={tt (B.2) 

0, h s = 0, 

and let S N denote the set of all distinct lists S with \S\ = N. Within S N the lists S 
and S' are ordered by letting j = min(z : Hi ^ n'^ and then setting S > S' if rij > n'-. 
Finally if S £ 5^ and 5" £ «S Ar ' with N > N' then S > S'. It is convenient to denote by 
5 + 1 the lowest list above S, and by S U S" the list obtained by concatenating S and 
5" and then ordering as above. 
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Now define 

H(S;f(x)) = (-if 1 n/ (ni) (^)> (B-3) 

i=l 

and for the empty list H([]\ f(x)) = 1. We need the following lemma, which is proved 
in 



Appendix B.l 



Lemma 

(i) If (H(S; Pt(x))) is finite as x — > for all S < S then (H(S; Pc{x))) is finite as 
x — > for all S < S and Ei ^ 0. 



(ii) If the conditions of part (i) apply and, as x — > 0, \H(S + l,Pr(x))) diverges as 
x~ 7 ,7 > ; £/ien (H(S + l,Pc(x))} diverges as x~ 7 . 

Differentiating times gives 



Wc 1 j ~ (l-P c (x)) 2 + l j 2^ (l-P c (x))%+i lB ' 4j 

where C(S') is a combinatorial coefficient. It is straightforward to check for any 5 
that (H(S; Pt{x))) ^ a is analytic for e 2 7^ 0, and that (H(S; Pt{x))) B is analytic when 
e 2 < 0. When e 2 = 

P,(x))U =0 = c 5 £ 2 ' 5 '-^ (1 + 0(r 2 )) (B.5) 

from which (H(S; P^{x))) B is divergent for S = [[a/2]], and with smaller degree for 
[[a/2] — 1, 1] if 2k < a < 2k + 1, k 6 Z, but always convergent for any inferior S 1 . The 
results given in section H~2l then follow from noting that P*(x) < Pc{x) < P(t(0) < 1 and 
using the lemma. 

Appendix B.l. Proof of lemma 
To prove the lemma note that 

H(S;f + g)= Yl H(S';f)H(S";g) (B.6) 

S'US"=S 

and differentiate ([HD n times to get 



where 



l) n P^ J (x) = (1 - x)F^(x) + nF'c \x) (B.7) 

SG<S n v v 

£ f^) /lS "^;^ 1 (x))P(S";P Tl (x)). 



.-.(«) 



s'us"=s 



(B.8) 
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It is then straightforward to generalise this formula to 



H(S,Pc(x))=K+(P c (x)) hs C(S,S')(-^ c{x)b 



s' 

X 



x )b- 

y) H(S";P Cl (x))H(S'";P Tl (x)), 



S'<S 
he 



s"us"'=s' 

(B.9) 

where the leading terms are written out explicitly and TZ contains contributions 
depending only on lists inferior to iS' 5 '. Every term on the right hand side is positive so 
it can be bounded above by using Pc(x) < -Pjt(O) and then the expectation value taken; 
moving the S" = S term to the left hand side gives 



(H(S,P c (0)))^(l-(^^" hS 



< 



n+(P t (o)) hs E c & s '^ 



S"uS"'=S' 

s'Vs 



s'<s 



6_ 



X 



(B.10) 

Part (i) is true for S = [1] so the lemma then follows immediately by induction on 5*. 
To prove part (ii) use part (i) to isolate the potentially divergent terms in (IB. 91) leaving 

H{B,P c (x))= { j0^f ) hS (h(S;P Ci (x)) + (^j hS H(S;P Tl (x)) 

+ finite terms. (B-H) 

For small enough x, 

«<(^)<^VC (B, 2 ) 
and part (ii) follows upon taking expectation values. 

Appendix B.2. E\ = 0, e 2 > 

We will show that 

_ / H(S,P c (x)) \ 
Fs -\(l-P cW )^/^ (R13) 

is finite at x = 0, which together with (1B.4I) gives the result. Using (loTT) and (1C.9I) gives 
F s <(n h c s+1 H(S,P c (x))) . (B.14) 
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where nc is the location of the first infinite tooth of C. Applying (1B.9P iteratively we 
find that the right hand side is bounded above by terms of the form 

(n£>^(#(^P T (z))>^. (B.15) 

The maximum value of K occurring is hg + 1 + 3>s where $s is the number of strings 
inferior to S. As remarked before (H(S', Pt(x))) ^ a is analytic and (^c)^a is trivially 
finite which completes the proof. 



Appendix C. Calculation of Pc(x) for some useful combs 

Let the comb C have T k = Ni,k < L and arbitrary Tl and Cl- Then following the 
method of Appendix A of [14J we find 

/ (1 - A)(P" 162 (x) - P h T 2 (x)) \ 

where 

(1 - x)b. 

{p;n*)?K (a2) 

Setting e 2 = 0, e\ < 0, £ = 1, Tl = Noo and Cl = tl we find after some algebra that 
Pf(x) = Jf V) (l + A- L + 0(*)) ) (C.3) 

and, as x — > 0, 

^ - i±N • (C.4) 

1 ~~ l e l| 

Repeating the exercise but with Cl = * yields a similar result. 

If instead we set e 2 > 0, e\ < 0, £ = 1, Tl = Noo and Cl = (j we find 

Pg*(x) = If»(x) (l + ^'^l^ (! + °(^))) ( C - 5 ) 
and, as x — > 0, 

1 + l , ei| ~ e2 . (C.6) 
1 - |ei| - e 2 

Again, repeating the exercise but with Cl = * yields a similar result. 
With e 2 > 0, e x = 0, and C = {T k<L = N e , T k > L = N^} we find that 

P^(x) > P^(x) - j^, L>2, (C.7) 
(it is good enough to use P*{x) for k = 2); and for C = \T k ^ L — Ni,T L — N^}, x < x , 

P**{x) < P^(x) I 1 - AL _^ 1 _ . L-1 j (C. 



where A = (1 - x)(P fl Ue2 (x))- 2 and 5 is a positive constant depending oni , A and e 2 . 
Finally for e 2 > 0, e\ — 0, and C = {T k ^L = N\, Tl = Noo} we find that 

W - 1 - L + v-«)/w (a9) 



Biased random walks on random combs 



26 



Appendix D. Upper bound on P(x) when ei = 0, £2 > 

We start by writing 

H{x-n) = {D^x))lp^-\ , (D.l) 



»\1-Pc{x). „ 

where the measure f2 is defined in ( )92|) . Applying Jensen's inequality with this measure 
to (J6]) results in a recursion formula of the same form as discussed in Appendix C and 
it is easy to verify that (P Ck (x)}- > (P Ck (x)) ^ to give 

G c (x;n) \ I G c {x;n) 



1-P c (x)/- ~ \1-P c (x 

3 exp (-n (1 - Pc(x)} ) 

> ^ ZL en 2) 

" b + (l-x)^ (l-Pcix))^ ' 1 • 1 

where in the last line we have again used Jensen's inequality when averaging over the 
ensemble. Applying this result to (ID. 1[) . summing over n, and using (1241) we obtain the 
inequality 

->(D\ T dx)) 2. (D.3) 

x~\ |T " 11 ^ (1 - P c ( x ))l 

Using the value for (Di(x)) given in tableland rearranging gives the upper bound on 
P{x) quoted in 15.2.31 



References 

[1] G. H. Weiss and S. Havlin, Some properties of a random walk on a comb structure, Physica 134A 
(1986) 474-484 

[2] S. Revathi, V. Balakrishnan, S. Lakshmibala and K. P. N. Murthy, Validity of the mean-field 

approximation for diffusion on a random comb, Phys. Rev. E 54 (1996) 2298-2302 
[3] D. ben-Avraham and S. Havlin, Diffusion and reactions in fractals and disordered systems, 

Cambridge University Press, Cambridge (2000) 
[4] S. Havlin, J. E. Kiefer and G. H. Weiss, Anomalous diffusion on a random comblike structure, 

Phys. Rev. A 36 (1987) 1403-1408 
[5] J. Ambj0rn, B. Durhuus and T. Jonsson, Quantum geometry: a statistical field theory approach, 

Cambridge University Press, Cambridge (1997) 
[6] J. Ambj0rn and Y. Watabiki, Scaling in quantum gravity, Nucl. Phys. B445 (1995) 129-144, 

[hepTh/9501049| 

[7] J. Ambj0rn, J. Jurkiewicz and R. Loll, Spectral dimension of the universe, Phys. Rev. Lett. 95 

(2005) 171301, |hei>th/0505113| 
[8] C. Aslangul, P. Chvosta and N. Pottier, Analytic study of a model of diffusion on a random 

comblike structure, Physica A 203 (1994) 533-565 
[9] S. Havlin, A. Bunde, H. E. Stanley and D. Movsholvitz, Diffusion on percolation clusters with a 

bias in topological space: non-universal behaviour, J. Phys. A 19 (1986) L693-L698 
[10] V. Balakrishnan and C. Van den Broeck, Transport properties on a random comb, Physica A 217 

(1995) 1-21 

[11] S. Havlin, A. Bunde, Y. Glaser and H. E. Stanley, Diffusion with a topological bias on random 
structures with a power-law distribution of dangling ends, Phys. Rev. A 34 (1986) 3492-3495 



Biased random walks on random combs 



27 



[12] N. Pottier, Diffusion on random comblike structures: field-induced trapping effects, Physica A 216 
(1995) 1-19 

[13] N. Pottier, Analytic study of a model of biased diffusion on a random comblike structure, Physica 
A 208 (1994) 91-123 

[14] B. Durhuus, T. Jonsson and J. F. Wheater, Random walks on combs, J. Phys. A 39 (2006) 1009- 
1038, |hep-th/0509i9ll 

[15] B. Durhuus, T. Jonsson and J. F. Wheater, The spectral dimension of generic trees, 
[math-ph/0607020 

[16] B. Durhuus, T. Jonsson and J. F. Wheater, On the spectral dimension of generic trees, DMTCS 

proc. AG (2006), 183-192 
[17] T.M. Elliott, Oxford University D.Phil Thesis, in preparation. 

[18] W. Feller, An introduction to probability theory and its applications, Vol.2, Wiley, London (1968) 
[19] A. Grigoryan and T. Coulhon, Pointwise estimates for transition probabilities of random walks 

in infinite graphs, in: Trends in mathematics: Fractals in Graz 2001, Ed. P. Grabner and W. 

Woess. Birkhaueser (2002) 



